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Abstract 

We investigate the effects of strain on a crystal surface close to the bulk melting 
temperature Tm, where surface melting usually sets in. Strain lowers the bulk 
melting point, so that at a fixed temperature below but close to Tm the thickness 
of the quasi-liquid film is expected to grow with strain, irrespective of sign. In 
addition, a strain-induced solid surface free energy increase/decrease takes place, 
favoring/disfavoring surface melting depending on the sign of strain relative to 
surface stress. In the latter case one can produce a strain-induced prewetting 
transition, where for increasing temperature the liquid film suddenly jumps from 
zero to a finite thickness. This phenomenology is illustrated by a realistic molecular 
dynamics simulation of strained Al(llO). 

KEYWORDS: Surface thermodynamics, Surface melting, Surface stress, Prewet- 
ting, Aluminum, Molecular dynamics 
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1 Melting of a Strained Surface: Phenomeno- 
logical Theory 



We consider the fate of the surface of a strained sohd just below the bulk 
melting temperature, and focus in particular on surface melting, that is on the 
possibility that a microscopically thin surface film could melt before the bulk. 
The phenomenological surface free energy variation (per unit area) produced 
by melting the solid surface, subject to a small bulk strain e (assumed parallel 
to the surface), to a thickness / of liquid can be written as 

AG = IpX (l - + ilsL + ILV - isv) + Vil) - ^Y'le' - a^'h . (1) 

Here A is the enthalpy of fusion per unit mass and p is the liquid density; 
thus the first, positive term represents an increase of free energy (T < Tm) 
required for melting a surface unstrained solid film to become a liquid film 
of thickness /. The second, A700 = {'Jsl + Ilv — Isv)^ is the free energy 
variation caused by replacing, at zero strain, the solid-vapor interface with 
the solid-liquid and the hquid-vapor pair of interfaces. When A700 < 
ordinary, strain-free surface meltingfll] takes place, with / > close to Tm- 
The third term V{1) represent effective interaction, usually repulsive, between 
the solid-liquid and the liquid-vapor interfaces. Its typical behaviour for 
small / is V{1) = \A\e^~'^''^^^ (A is a parameter comparable with A700, ^ is 
the correlation length in the liquid 0), while a crossover to the asymptotic 
form V{1) = H/P appears at larger Z's. The fourth term represents the 
decrease of the elastic energy stored in the strained solid, {l/2)Y'le'^ , since 
the molten film is free to expand or contract along the surface normal {Y' is 
proportional and close to the Young modulus y)0. Finally a'-'^^e is a surface 
free energy change - in principle affecting both solid-vapor and solid-liquid 
interfaces - caused by a generally nonzero surface stress o"'-'*^ . 

Eq.(|l|) can be rewritten as 

AG(/) = p\ {^^^M_T\j I - (IA700I + <y^'^e) + V{1) , (2) 

^ The exact value of Y' depends on how the strain is apphed. If the surface is strained 
in the y direction [syy = e) and the x direction is free (uxx = 0) in the absence of shear 
(Sxy = 0), then Oyy = Ysyy and Y' = {p/paoiid)Y- The density ratio enters because /, the 
thickness of the molten film, is larger than the thickness of the original solid film. 
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where T^, is defined as 

As the temperature reaches T^, the thickness of the molten layer diverges: 
T^{e) is the melting temperature of the solid under strain e. Actually the 
strain enhances the free energy of the solid phase without altering that of 
the liquid phase. Thus the latter is favoured and < T^- We note that 
AG(0) = but limj_o+ AG(/) = G(0+) = -(IA700I + a^'h) + ^(0+) which 
does not generally vanish. In fact equations (|I]) and (0) cease to hold for / 
comparable with a monolayer or less. 

The predicted strain-temperature phase diagram for realistic parameters 
mimicking aluminium, and obtained minimizing the AG(/) with respect to I 
is shown in Fig. |1]. At T = T^{e) (solid line) there is melting of the strained 
bulk. The quadratic decrease of with increasing strain is visible. Sur- 
face melting appears below when, for increasing T, AG(/) first develops 
a minimum for a finite value of The nominal onset temperature for sur- 
face melting - the temperature where the solid surface is first wetted by 
an infinitesimal liquid film - is attained at T = where d{AG)/dl = 
at / = 0. This wetting temperature (dot-dashed line) is, similarly to bulk 
melting, quadratically depressed by strain. 

Between these two temperatures, and T = T^, the quasi-liquid film 
thickness grows from zero to infinity. The divergence at is power law 
/ ~ (T^ - T)-i/3 H > 0, but only logarithmic / ~ \logiT;^ - T)\ for 
H = 0. A crossover between the two regimes appears when the derivative of 
H/l"^ with respect to / becomes comparable with that of | A7oo|e~^'/'^: in Al 
H ~ 0.56 X 10~^^ Joule[Q] and this crossover should take place at about 0.5 
K below the melting point, with / ~ 14A. 

A continuous growth of the liquid film thickness is typical of regular, 
complete wetting of the solid substrate by its own melt. For other surfaces 
or materials wetting may not occur at all before bulk melting (T^ = T^), and 
we have the so-called surface non-melting[|, ||, A discontinuous growth, 
with a jump in the liquid thickness, could theoretically take place in case of 
a so-called prewetting |^ transition. The strain-free surfaces that have been 
studied so far, experimentally as well as theoretically, were found to exhibit 
mostly complete surface melting, or surface nonmelting. An intermediate 
type of behavior known as incomplete surface melting, where the liquid film 
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Figure 1: The phase diagram of the Al(llO) interface. The circles show 
the melting temperatures of the strained solid obtained by the molecular 
dynamics simulations. The prewetting line has been obtained by assuming 
V{1) = I A7oo|e~^'/'^ in equation (|l]). Note the large strains we are considering. 
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thickness levels off to a finite value below T = T^, disco ntinuously jumping 
to infinity at T = T^, has also been described § and can be seen as a rather 
special case of prewetting. However no case of regular prewetting, with a 
finite jump of the liquid film thickness, has so far been described in surface 
melting. 

As it turns out, our simple model also predicts that surfaces that exhibit 
complete surface melting could be caused to develop a prewetting transition 
by means of external strain. That is due to the inevitable presence of nonzero 
surface stress at the solid surface, or more correctly at all interfaces involving 
the solid. The ere term provokes a shift, linear to first order in the strain, 
of the interface free energy balance A7oo. Depending on that shifted value 
the minimum of AG{1) can be either negative with respect to the crystalline 
surface - thus supporting a stable liquid film - or positive, in which case 
the liquid film is metastable, supporting a stable dry solid surface. In that 
case there is a whole range of strain values where AG{1) changes its sign 
from positive to negative for increasing temperature, generating a prewetting 
transition phase line, where I jumps from zero to a finite value. For even 
larger negative strain magnitudes the surface behavior is eventually predicted 
to become nonmelting. 



2 Simulations of Strained Surface Melting 

To verify the above simple theory, we simulated the thermal behavior of 
Al(llO) close to the melting point, in presence of unidirectional in-plane 
strain. Molecular dynamics simulations of Al(llO) were done in the slab 



geometry, both flat and bent, ||T0[ and the Ercolessi-Adams many body po- 
tential for Al, whose bulk melting temperature is 943 K (to be compared with 
an experimental of 933 K) was used. Strain was introduced by expanding 
the flat simulation box along the [0 1] direction, or alternatively by bending 
the simulated slab, while keeping it fixed in the orthogonal direction. 

The simulations, conducted with standard canonical methodsconfirmed 
first of all the decrease of T = T^[e) with strain. Fig. ^ shows the large 
increase of the melted film occurring at fixed temperature for symmetrically 
positive and negative strains. 

Subsequently, to investigate the possible presence of a prewetting tran- 



sition, we carried out iTTIl a series of microcanonical simulations, where at 



each given strain the internal energy was increased stepwise. (Fixing the en- 
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Figure 2: Single snapshots of the simulations of Al(llO): side view of of three 
samples at the same temperature. Both negative and positive strain increase 
the molten film thickness, reflecting lowering of the melting temperature by 
strain. Samples size: 14 x 20 x 16 cells. The atom stacking in the 20-atom 
rows orthogonal to the picture clearly distinguishes the liquid by the solid 
phases. 
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Table 1: The decrease of the effective melting temperature with strain 
The third column is obtained by formula (0), with Y'/{2pX) = 26.6 for the 
Ercolessi-Adams potential. 



ergy in place of the temperature is useful to reduce the large fluctuations of 
constant temperature simulations close to a first order phase transition). At 
each step the temperature as well as various structural correlation functions 
were monitored along 2.12 nanosecond long runs. The effective melting tem- 
perature is extracted, by determining the vertical asymptote in the plot 
of temperature versus molten film thickness. Data reported in Table ^ show 
a good agreement with our formula (Q), the discrepancy at the larger posi- 
tive strains being due to anharmonic effects. In this kind of simulation the 
possible occurrence of prewetting will show up, as with any other first order 
transition, by hysteresis, as well as by the occurrence of an inhomogeneous 
two-phase coexistence of dry and wet portions of the surface. Because the 
surface stress of Al(llO) is positive (here about 0.052 eV/A^ at a temperature 
of 900 K, but ab initio calculations of Needs [jl2| at T = provide a value of 
0.115eV/A^), it is expected that prewetting could appear in an important 
way at negative strain, i.e., under compression. 

Fig. 1^ shows the effective liquid thickness against temperature as obtained 
in each simulation, for strains of -3%, 0, and +3%. The expected overheating 
of the solid surface is indeed observed for strains -3%, 0, but not for +3%, in 



agreement with the snapshots. We conclude, as detailed in Ref. ||TT[, that the 
Ercolessi-Adams model potential predicts a prewetting transition in the sur- 
face melting of compressed and of strain free Al(llO). The prewetting jump 
in the liquid layer thickness disappears form the huge positive strain of +3%. 
The agreement with Fig. |I| is qualitatively correct, although the prewetting 
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Figure 3: Thickness of the molten layer vs Temperature. Continuous and 
dashed lines are the predictions, without considering the prewetting, for un- 
strained and strain e — ±0.03 slabs respectively. 



8 



region turns out to be larger than expected, unexpectedly including zero 
strain. 



3 Discussion and Conclusions 

Bulk strain depresses the melting point of a solid, whose surfaces will tend, at 
a given fixed temperature, to melt more readily than in the absence of strain. 
Since a tiny strain of order 0.0006 suffices in Al to lower by one hundredth 



of a degree - the typical accuracy of a surface melting experiment [p!^] - this 
suggests that random strains could represent an important source of uncer- 
tainty in these experiments. A combination of strained regions, resulting e.g. 
from surface treatments, could lead to patchy behaviour close to melting, 
where portions of the surface would melt and others would not. 

Another consequence of strain-induced surface melting is the related pos- 
sibility to give rise to potentially interesting side effects. If strain were in- 
troduced through a bulk longitudinal wave for example of sufficiently low 
frequency u and wavevector /c, one could expect a corresponding surface liq- 
uid thickness modulation of frequency 2uj and wavevector 2k. There would 
also be associated frictional damping effects worth investigating. 

A newer effect of strain is the possibility to cause a prewetting transition 
in surface melting. That is of interest, in view of the fact that prewetting has 
never so far emerged in surface melting. Although it might prove difficult 
for the crystal to sustain large strains particularly so very close to the melt- 
ing point, nevertheless the entirely theoretical possibility of strain-induced 
prewetting seems worth addressing experimentally, as the stress needed might 
in fact be quite small. In the case of Al(llO) used here as a test case, our sim- 
ulations indicated prewetting already at zero strain. While that may be an 
artifact of the potential used - existing experiments have shown no evidence 
of prewetting on A1(11O)0 - it seems possible that prewetting could appear 
with a relatively modest compression, as suggested by Fig. |l[ Prewetting is 
further discussed in Ref. [0. 
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